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Abstract We present a first analysis of a deep X-ray 
spectrum of the isolated neutron star RBS1223 obtained 
with XMM-Newton. Spectral data from four new mon- 
itoring observations in 2005/2006 were combined with 
archival observations obtained in 2003 and 2004 to form a 
spin-phase averaged spectrum containing 290000 EPIC- 
pn photons. This spectrum shows higher complexity than 
its predecessors, and can be parameterised with two Gaus- 
sian absorption lines superimposed on a blackbody. The 
line centers, E2 — 2Ei, could be regarded as supporting 
the cyclotron interpretation of the absorption features 
in a field B ~ 4 x 10 13 G. The flux ratio of those lines 
does not support this interpretation. Hence, either fea- 
ture might be of truly atomic origin. 

Keywords Stars: individual (RBS1223) • stars: 
neutron • stars: magnetic fields ■ X-rays: individuals 
(RBS1223) 
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1 Introduction 

RBS1223 is one out of three XDINs (X-ray dim iso- 
lated neutron stars) found in the ROSAT Bright Sur- 
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vey (RBS), an optical identification program of the more 
than 2000 X-ray sources detected in the RASS with a 
count rate > 0.2 s -1 at high galactic latitude pp. The ini- 
tial discovery [2] was based on the soft X-ray spectrum 
and the steep SED, fx/ f opt > 10 4 - Follow-up Chan- 
dra observations revealed a periodically modulated sig- 
nal 3 , and HST-observations uncovered a candidate op- 
tical counterpart at msoccD = 28.56 ± 0.13 mag 0]. Ini- 
tial observations with XMM-Newton showed deviations 
from a Planckian energy distribution at energies below 
500 eV which could be described with a Gaussian ab- 
sorption line and interpreted as a cyclotron absorption 
line in a field of (2 — 6) x 10 13 G (SJ. The large number 
of photons collected with XMM-Newton also uncovered 
the true spin period of 10.31 s. In the light curves 
in different energy bands were modeled in terms of a 
photospheric cap model. The temperature structure was 
found to be roughly compatible with the crustal field 
models of The implied field structure deviates from 
a simple centered dipole model. Spin-period changes ap- 
peared likely at that time thus stimulating monitoring 
observations with Chandra and XMM-Newton. The spin 
history could be unequivocally fixed by a series of Chan- 
dra pointings |S] and revealed a spin-down of the star at 
a rate P = 1.120 x 10 _13 ss _1 . Under the assumption 
that this is due to magnetic dipole torques, a charac- 
teristic age of 1.5 Myr and a magnetic field strength of 
3.4 x 10 13 G were inferred. Here we describe initial results 
of recent monitoring observations with XMM-Newton 
which, due to the large number of accumulated photons, 
besides accurate timing also allows spectroscopy with un- 
precedented accuracy. 



2 XMM-Newton observations of RBS 1223 

XMM-Newton observed RBS1223 several times between 
2001 and 2006. The relevant pointings for this contri- 
bution are listed in Tab.Q] All but the first, which was 
put in small window mode (SW), were obtained in full 
frame mode through the thin filter. All observations were 
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Table 1 XMM-Newton observations of RBS1223 analysed 
in this contribution 



Rev. 


AO 


Day 


Exp. 


Mode 


Countrate 








M 




Is" 1 ] 


377 


1 


2001 Dec 31 


18.6 


SW 


2.49 a 


561 


2 


2003 Jan 1 


27.0 


FF 


2.53 


743 


3 


2003 Dec 30 


30.2 


FF 


2.57 


1015 


4 


2005 Jun 25 


14.9 


FF 


2.56 


1016 


4 


2005 Jun 26 


12.9 


FF 


2.57 


1025 


4 


2005 Jul 15 


12.9 


FF 


2.56 


1115 


4 


2006 Jan 10 


14.9 


FF 


2.56 



- Exposure time is given for EPIC-pn. 

- Mean countrates are given for the whole XMM-Newton 
energy range (0.15 - lOkeV) 

- a countrate from [5] 



reduced in a homogeneous manner with a recent ver- 
sion of the SAS (V6.5). The monitoring observations set 
constraints on the evolution of different parameters of 
the star. The overall brightness did not change between 
2001 and 2006 and at all occasions the star displayed 
its double-humped light curves discussed in detail by 
[7] . Also the amplitude of the light curve remained un- 
changed. Hence, the only property which was shown to 
evolve so far is the spin period [Sj. 



3 Spectral analysis of RBS1223 

The constancy of the star allows a joint spectral anal- 
ysis of all available observations. In order to minimise 
possible remaining calibration uncertainties for different 
camera modes we used only those observations with iden- 
tical camera settings (six observations between revolu- 
tions 561 - 1115). This selection left us with an effective 
exposure time of 113 ks with almost 290000 photons for 
spectral analysis. Here we concentrate on the mean spec- 
trum of all observations. Phase-resolved spectroscopy is 
underway and will be presented elsewhere. 

An initial attempt to analyse the spectrum was made 
with a spectrum extracted from a merged photon event 
table of the six observations and with common response 
and effective area files. Such an approach does not take 
into account possible time-dependencies of the detector 
files. We therefore extracted a spectrum and computed 
the response matrix and effective area file for each obser- 
vation separately. Spectral modeling was performed with 
XSPEC V12 following two different approaches. Firstly 
the spectral parameters were determined for each spec- 
trum individually but forcing the column density of cold 
interstellar absorbing matter, Nh, to be the same for all 
the observations (six data groups). This experiment re- 
vealed the spectral parameters of all observations to be 
the same within the errors. We then followed the second 
approach fitting all six spectra jointly with one common 
model which reveals one set of spectral parameters in- 



stead of six. The results of those latter experiments are 
listed in Tab. 

Spectral models were chosen with increasing com- 
plexity until a satisfactory fit was achieved. The fit as- 
suming only a blackbody emission spectrum absorbed 
by neutral interstellar matter which gives xt — 4 (for 
1247 d.o.f.) can be ruled out completely (model 1 in 
Tab. [21 Fig. [TJi). This mismatch was found already by 
[S] when fitting data from revolutions 377 & 561. They 
successfully applied a model with a Gaussian absorp- 
tion line superposed on the absorbed blackbody. With 
such a model the current fit is improved but it does no 
longer adequately represent the shape of the combined 
spectrum (model 2 in Tab.0 Fig.QJj)- Compared to [H] 
the fit with the single Gaussian gives a higher line en- 
ergy, Ei = 390 eV vs. 300 eV, and smaller line width, 
cti = 60 eV vs. 100 eV. Part of this change is due to 
the improved signal/noise of the spectrum which allows 
a free fit of o\ instead of fixing it at 100 eV as in 0. 
The updated calibration of the detector might also help- 
ful. After inclusion of a second Gaussian a successful fit 
showing no systematic residuals was achieved (xt = 1-11 
for 1244 degrees of freedom, model 3, Fig.QJ:). 

Since the second Gaussian could be the higher har- 
monic of a fundamental at lower energies (py proposed 
an interpretation as proton cyclotron line in a field of 
(2 — 6) x 10 13 G) we made some attempts to fit the two 
Gaussians with related parameters. In our model 3, both 
Gaussians were fitted independently, they were centered 
at line energies E2 — 2 x Ei . Hence, if the parameters 
are forced to obey this relation, the achieved fit becomes 
equally good and the other parameters are found almost 
unchanged (model 4) . The picture changes slightly if the 
width of the two lines is pre-set to the same, adjustable 
value (model 5). This broadens the width of the princi- 
pal first line slightly from 0.15 keV to 0.17 keV and shifts 
the line energy of the second Gaussian from 0.46 keV to 
0.73 keV at the same quality of the fit (xl = 1-H)- Nev- 
ertheless the flux in the second Gaussian is reduced by 
50% in this fit. The last model listed in Tab. num- 
ber 6) uses a different relation between the line energies, 
E2 = 1.5 X Ei, and was similarly successful with almost 
unchanged parameters compared to the free fit. 

For models 1 and 2 no formal parameter errors were 
determined because those models did not adequately de- 
scribe the spectrum. The formal errors of the parame- 
ters of fit 3 are: AN H = 0.2 x 10 20 cm- 2 , AkT = 2eV, 
AEi = 0.02 keV, Aa x = 0.02 keV, AF 1 = 4 x 10~ 4 keV 
cur 2 s"\ AE 2 = 0.03 keV, Aa 2 = 0.03 keV, and AF 2 = 
4 x 10~ 4 kcV cm~ 2 s _1 . The above experiments have 
shown cross-talk between the parameters, hence the true 
uncertainties are larger than the mentioned statistical 
ones. For example, models 3-6 are indistinguishable 
on statistical grounds but some of their parameters are 
rather different. The main problem consists in a proper 
location of the continuum. The primary Gaussian is so 
broad, that the low-energy end of the spectrum is still af- 
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Rbs 1323 XMM-Newton EPIC PN Rbs 1223 XMM-Newton EPIC PN Ebs 1223 XMM-Newton EPIC PN 



Model: phabs ( bbodyrad ) Model: phabs ( gaussian + bbodyrad ) Model: phabs ( gaussian + gaussian 4 bbodyrad ) 




0.2 0.5 1 0.2 0.5 1 0.2 0.5 1 



Energy (keV) Energy (keV) Energy (keV) 

Fig. 1 Results of spectral fitting ordered according to increasing complexity of the model. From left to right the models 
shown are an absorbed blackbody, a blackbody plus Gaussian absorption line, and a blackbody with two Gaussian absorption 
lines. The spectral parameters are listed in Tab.|21 models 1-3 



Table 2 Spectral modeling of RBS1223. E and a denote the central energy and width of the two Gaussians (subscript 1 
and 2), while F denotes the line flux 



# 


N H 


kT bh 


Ei 




Fi 


E 2 


02 


F 2 


red. x 2 


#d.o.f 




10 20 cm" 2 


eV 


keV 


keV 


keV cm -2 s -1 


keV 


keV 


keV cm~ 2 s _1 






1 


3.3 


100 














3.96 


1250 


2 


3.7 


93 


0.39 


0.06 


-6 x 10 -4 








2.05 


1247 


3 


1.8 


102 


0.23 


0.15 


-35 x 10" 4 


0.46 


0.26 


-7 x 10" 4 


1.11 


1244 


4 


2.0 


104 


0.28 


0.13 


-21 x 10~ 4 


= 2 x E 1 


0.24 


-7 x 10" 4 


1.12 


1245 


5 


1.2 


100 


0.20 


0.17 


-27 x 10~ 4 


0.73 


= O"! 


-3 x 10" 4 


1.11 


1245 


6 


1.8 


102 


0.28 


0.15 


-35 x 10~ 4 


= 1.5 x Ei 


0.26 


-7 x 10 -4 


1.11 


1245 



fected by this feature. This allows the flux in this features 
to be traded against the amount of interstellar absorp- 
tion and the temperature of the blackbody. At the high- 
energy end of the spectrum photons are detected up to 
2 keV which in principle gives a better leverage for con- 
straining the blackbody temperature. Our fitted model 
assumes that these photons belong to the Wien tail of the 
blackbody. Consequently, the newly determined black- 
body temperature is higher than found before from a 
spectrum which missed those high-energy photons and 
was affected by the structure which is now described by 
a second Gaussian in its presumed Wien-tail. 



4 Results and discussion 

The joint analysis of six observations of RBS1223 with 
XMM-Newton has led to a slight revision of the spec- 
tral parameters. The spectrum can be described with a 
blackbody and two Gaussian absorption lines superim- 
posed. The blackbody temperature is significantly higher 
than previously found thanks to the detection of X-ray 



photons up to 2keV. Also the parameters of the pri- 
mary Gaussian absorption line are changed: it appears 
significantly narrower and more blueshifted than before. 
Whether we have seen the true continuum at any energy 
remains an open question. 

Similarly unclear is the nature of the features which 
were parameterised with two Gaussian absorption lines. 
They have a flux ratio of about 5:1, their equivalent 
widths of EWi - 200 eV (for E 1 = 0.23 keV) and EW 2 ~ 
180 eV (for E 2 = 0.46 keV) are rather large. Their rela- 
tive line energy and the rather large width of the lines 
are supportive of the proton cyclotron interpretation in 
a field of few times 10 13 G (the line center of the pri- 
mary and hence the derived magnetic field strength is 
not well contrained, in model 5 the best fit is achieved for 
E\ = 0.2 keV, i.e. it converges to the lower bound of the 
pre-set parameter range), consistent with the observed 
spin down rate of the star. The rather high flux ratio of 
the two lines, nevertheless, makes an interpretation of the 
two lines as being harmonics of each other unlikely, since 
the oscillator strength of the first harmonic is smaller 
than that of the fundamental by a factor ~ E/(mc 2 ), 
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which becomes very small for proton masses (2). An al- 
ternative could be to associate the two lines with differ- 
ent fields from e.g. the two polar caps. Since the field was 
shown to be somehow non-dipolar or off-centered at least 
jZj, different field strength could be encountered at the 
two caps. It is difficult to test this hypothesis on the basis 
of the current data and it needs to be seen if the results of 
the phase-resolved spectroscopy reveal any helpful indi- 
cation in this respect. However, if the two lines would be 
the fundamentals of two different cyclotron systems from 
two different regions, the infered field strengths would be 
B x = 4.2 x 10 13 G and B 2 = 8.5 x 1 13 G according to 
B(10 13 G) = I ig(l + ^)2E^(keV) and assumed z = 0.18 
(see [7j), Ei = 0.23 keV and E 2 = 0.46 keV. Hence, one 
field would be slightly below and one above the criti- 
cal quantum field, where vacuum polarization becomes 
important in shaping the lines (see J^l j an d the dis- 
cussion in EI])- Since the second Gaussian line (higher 
field) is rather broader than the primary Gaussian, this 
scenario seems to be unlikely. 

A possibility which needs to be explored in detail is 
a blend of magnetically shifted atomic transitions with 
proton cyclotron resonances. As the model computations 
of |12| indicate, the amount of ionisation in the assumed 
hydrogen atmosphere will play a crucial role. An as accu- 
rate as possible temperature determination is thus very 
important for further progress. Due to its photometric 
stability on the one hand and the pronounced spin vari- 
ability on the other hand, RBSI223 is a very promising 
target to gain further insight from even deeper X-ray 
spectroscopy. 
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